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In this paper, we study the effects of instrumental systematics on the reconstruction of the deflection angle 
power spectrum from weak lensing of Cosmic Microwave Background (CMB) temperature and polarization 
observations. We consider 7 types of effects which are related to known instrumental systematics: calibration, 
rotation, pointing, spin-flip, monopole leakage, dipole leakage and quadrupole leakage. These effects can be 
. . , characterized by 1 1 distortion fields. Each of these systematic effects can mimic the effective projected matter 

power spectrum and hence contaminate the lensing reconstruction. To demonstrate the effect of these instru- 
' mental systematics on CMB lensing measurements, we consider two types of experiments, one with a detector 

' noise level for polarization of 9.6 //K-arcmin and FWHM of 8.0', typical of upcoming ground and balloon-based 

CMB experiments, and a CMBPol-like instrument with a detector noise level for polarization of 2.0 yuK-arcmin 
' and FWHM of 4.0', typical of future space-based CMB experiments. For each systematics field, we consider 

I various choices of coherence scale o-j, starting from a, = 10' to the maximum of a, = 120'. Among all 

, the 1 1 systematic parameters, rotation o), and monopole leakage ja^Tb place the most stringent requirements, 

' while quadrupole leakage q, pointing error pa and p^, and calibration a are among the least demanding. The 

requirements from lensing extraction are about 1-2 orders of magnitude less stringent than the requirements 
I I , to measure the primordial B-modes with inflationary energy scale of 1.0 x 10'* GeV. On the other hand the 

■ requirements for lensing reconstruction are comparable or even more stringent for some systematic parameters 

' than the requirements to detect primordial B-modes with inflationary scale Ej = 3.0 x 10'* GeV. 
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I. INTRODUCTION 



Observations of the temperature anisotropies of the Cosmic Microwave Background (CMB) have been a major tool to con- 
^ strain cosmological parameters. The polarization data of the CMB can help us to extract additional information beyond the 
temperature information Oli^. The next generation of CMB observations will focus on the precise measurement of polarization 
00 of the CMB, especially the so called B-mode polarization, which is at least two orders of magnitude smaller than the tempera- 
04 ture anisotropy signal. In contrast to E-mode polarization, which can be generated by scalar or tensor perturbations in the early 
', universe, the primordial B-modes are generated only by tensor perturbations IHQl- However even in the absence of primordial 
T— I ■ B-modes, subjeguent^gravitational lensing by the large scale structure of the Universe converts E-mode polarization to B-mode 
' polarization 15|, |6|, |7j, [81] . Although the amplitude of the primordial gravitational wave signal is uncertain by many orders magni- 
tude and might not be detectable by the next generation of polarization experiments, the lensing B-mode signal is a guaranteed 
prediction of the current cosmological model. In addition, the B-mode lensing signal will help to break degeneracies between 
cosmological parameters |9, 10, 11, 12, 13, l^l- 

The weak gravitational lensing of CMB anisotropies provides a unique opportunity to map the matter distribution of the 
universe. The non-Gaussian higher order correlations in the CMB generated by the weak lensing can be used to reconstruct the 



mass distribution of the intervening large-scale structure. The principle is easy to understand. The CMB photons are remapped 
by gravitational lensing which introduces coiTelations between different angular moments. One way of extracting the lensing 
information is to use a quadratic combination of the CMB multipoles to define an estimator for the projected gravitational 
potential SEIBIilEIISIIllIl. 

Although the B-mode polarization observations are currently still noise dominated, the next generation of CMB polarization 
instruments has the sensitivity to make first detections, at least of the lensing induced B-mode signal. However, there are 
several challenges for CMB lensing detection, mainly coming from astrophysical foregrounds and instrumental systematics. It 
is important to estimate and control those spurious signals as well as possible when analyzing upcoming CMB data. These 
challenges will have to be overcome in order to probe the physics of the early universe through B-mode polarization or to infer 
the projected large scale matter distribution from E/B polarization. 
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Lensing studies can be considered as secondary science for an experiment devoted to B-mode detection. Impact of instru- 
mental systematics on the projected matter power spectrum is helpful to both instrument design and future data forecast. On 
the other hand, lensing induced B-mode polarization is a cosmological contaminant for the detection of primordial B-modes, a 
systemic study of instrumental systematics for lensing reconstruction may well be required of such an experiment to delense the 
observed CMB fields 1231 ^ \25, ,26, .27, ,28, ,2% ,30, ,3 Hi . 

In the literature, instrumental systematics have been discussed extensively iH^-SHIIsISIItIIIMII^- The main goal of 
this paper is to illustrates the effects of instrumental systematics and systematically study the impact on the mass reconstruction 
process for upcoming CMB experiments [24, 25, 26, 27, 28, 29, 30]. To calculate the effects of instrumental systematics on 
the projected matter power spectrum, we make use of the quadratic estimator to reconstruct the projected gravitational potential 
ifTsll . and classify and parameterize the systematics following ll32ll . We divide polarization contaminations into two categories: 
those which are associated with a transfer between the polarization state of the incoming radiation (from detection system), 
including calibration and rotation, spin-flip coupling and monopole leakage errors, and those which are associated with CMB 
anisotropy induced by the finite resolution or beam of the telescope. For the polarization transfer systematics, we only concern 
ourselves with polarization transfer in a single, perfectly known, direction on the sky. However, in reality, every experiment 
necessarily has finite resolution and this therefore is an additional class of contamination associated with the resolution or beam 
of the experiment. We refer the reader to ll32ll for a detailed discussion of the parametrization of the systematic errors we use 
in this paper We calculate the unlensed and lensed CMB power spectrum using CMBFAST ll40ll . In the calculation, we have 
assumed a flat ACDM cosmology with following cosmological parameter values: Q/, = 0.045, Qc = 0.23, //o = 70.5, n, = 
0.96,«, = 0.0, T = 0.08. 

This paper is organized as follows: In Sec.|ll]we review the basic lensing formalism, define our notation, and introduce the 
quadratic estimator of the lensing potential reconstruction which we later use to explore the effects of instrumental systemat- 
ics. In Sec. [nil we first consider a relatively simple case of the calibration systematics on the temperature estimator of lensing 
reconstruction. Then we consider EB estimator as an example, studying 7 types of instrumental systematic effects on lensing 
reconstruction, and compare to the systematic contamination of the B-mode power spectrum detection. We consider two instru- 
mental configurations, one with noise sensitivity for polarization of 9.6 //K-arcmin and FWHM of 8.0' (Expl from here on), 
and another CMBPol like instrument with noise sensitivity for polarization of 2.0 yuK-arcmin and FWHM of 4.0' (reference 
experiment from here on). In Sec. II VI and [Vl we describe our results and conclude with a discussion of the implications for 
experiments dedicated to measure primordial B-modes or the secondary lensing signal. We leave discussions of the other three 
lensing potential estimators (EE, TE, TB) to Appendix. 

II. LENSING FORMALISM 

Gravitational lensing deflects the path of CMB photons from the last scattering surface resulting in a remapping of the CMB 
temperature/polarization pattern on the sky. In this section, we review the basic lensing calculation for both temperature and 
polarization fields as the starting point of our discussion. We formulate CMB lensing using the flat-sky approximation ll22ll . The 
flat-sky approach simplifies the derivation by replacing summations over Wigner symbols of spherical harmonic moments by 
integrals involving mode coupling angles ||7[]. More details of CMB lensing can be found in the nice review paper ll4lll . 

At a certain position n on the sky, the observed CMB field X(n) is lensed from another direction in the primordial CMB sky 
X(n + d(n)) at z=1090. The remapping process can be described as 

f(n) = r(n + d(n)), (1) 
[Q + iUm = [Q±iU](A + dm, 

where 7'(n) (^(n)) represents the lensed (unlensed) temperature fluctuation field, Q{n) {Q{n)) and U(n) iU{n)) are lensed (un- 
lensed) polarization Stokes parameters, and d(n) is the deflection angle which is related to <p{fi), the lensing gravitational poten- 
tial, by d{n) - V0(n). Here and throughout this paper, we use boldface quantities to identify vectors, and X (X) stands for lensed 
(unlensed) temperature and polarization fields. The lensing potential 0(n) is given by 

m = -2 dr /\l ' 0(r,m), (2) 
Jo dAindAirQ) 

where dA is the comoving distance along the line of sight; ro is the comoving distance to the surface of last scattering, and O is 
gravitational potential. The lensing remapping process conserves the surface brightness distribution of the CMB, thus does not 
change the one-point statistics. 

The observed temperature and polarization fluctuations also include secondary effects, such as Sunyaev-Zel'dovich (SZ) 
effect and Integrated Sachs Wolfe (ISW) effect ll43ll . which come from the first order density or potential fluctuation and 
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thus also correlate with the lensing deflection angle. We denote these physical contamination to lensing reconstruction by X""^(n). 
We denote the noise component by X"{n). The total observed CMB anisotropy therefore includes the lensed primary signal, any 
secondary eff'ects, and noise, i.e. X'{n) - X{n) + X''"^(n) + X"{n). In the next section, we will introduce another contribution 
to X'(h), which comes from instrumental systematics X"''(n). We define the observed CMB field X"''^(h) = X'(h) + X"'\n). 
Here we write the secondary contribution as an independent component from the lensed CMB. However, in reality it is hard 
to separate X^'^'Xn) from X'{n) because secondaries are also lensed by gravitational potentials with deflection angles depending 
on their redshifts. In this paper, we simply drop the contribution effects X^'^''{n), as this topic is beyond the focus of this paper 
We refer the readers to | 22| l for a treatment of the secondary anisotropy as a physical contamination to the lensing potential 
reconstruction analysis'. 

It is convenient to work in Fourier space. If one considers a small enough patch of sky, spherical harmonic modes can be 
replaced by Fourier modes. Generalization from the flat-sky to the full sky is straightforward. The Fourier transform of the 
Taylor expended lensed CMB temperature and polarization field is 



), (3) 

21. 



til) = J flfnf(n)e-"* = r(l)- J ALr(l')L(l,l' 
[e ± iB] (1) ^ Jdh [Q{h) ± /J/(n)]e^2,V'e-'i" = [£(1) ± iB{l)] - J ^ 1^') ± iB(l')]Lp(l 1') , 

(^(1) = J dncf)(h)e-'^-^, (4) 



where 



L(l, 1') = ,^(1 - 1') [(1 - 1') ■ 1'] + i r ^0(1") X 0(1 - 1' - 1") d" ■ 1') [(!" + 1' - 1) . 1'] + . . . , (5) 

Lp(l, 1') = e^2.(w-w)0(i _ 1') [(1 _ 1') . 1'] + ^ ^(l _ l' _ l") (l" . l') [(l" + i' _ i) . i'] + . . . . 

2 J (2n)~ 

We can immediately see that lensing induces remapping of CMB fields by lensing potential gradients, hence in Fourier 
space lensing acts as a convolution which couples different harmonic modes. Fourier moments, power spectrum, bispectrum, 
trispectrum and so on of the CMB fields and the lensing potential can be defined in the usual manner: 

{x'(h)X'\l2)) = (27rf6u(h+h)C'^,^,Ah), 

{X'ili)X'\l2)X"'(h))^ = (27r)25D(li +I2 +l3)fi^X'X"(ll'l2'l3), 

{r{ii)x'j{h)x"'(h)x"""(h))^ = (2;r)25D(ii +I2 + I3 + i4)ri'^:^„;,,„(ii, 12,13,14), 

(6) 

where the angle brackets represent ensemble averages over realizations of the primordial CMB fields, the large-scale structure 
between observers and the last scattering surface, and the experimental noise. The connected part of the n-point function is 
denoted by the subscript c. The fields X, X', X" , X'" are among {^(l), £(1), Z?(l), 0(1)). The superscripts ;, j, k, m represent the 
unlensed field X, the lensed field X, the instrumental noise X", the instrumental systematics X'"'\ CMB secondary contribution 
X""'', the total signal X', or the observed field X°'"'. We note that the bispectrum and higher order odd-correlations vanish if one 
ignores the secondary effects. This is because odd moments contain sample averages over the odd primordial CMB fields which 
we assume to be Gaussian. 

We make the assumption that fluctuations in the large-scale structure between the observer and the last scattering surface are 
Gaussian and hence can be fully described by a power spectrum. We use the lensing potential power spectrum calculated from 
CAMBFAST. The instrumental noise X" is also assumed to be Gaussian. We note that primordial non-Gaussianity can possibly 
contribute percent level uncertainty to our analysis which is considered in |45!!|. 

We are now in the position to calculate any order (cross) correlation functions of CMB fields and the lensing potential in 
Fourier space. We will assume uniform Gaussian noise with the power spectrum Cf^" given by 

Cf^" = w^^e'^'-'i . (7) 



As a note, thermal SZ effect can in principal be separated from the primary fluctuations by its spectral dependence. For the kinetic SZ effect, it was claimed 
that by using a specially designed estimator |44|| . it is possible to separate it out from real lensing signal. However, some important secondary contributions 
such as the ISW effect cannot be separated easily and will lead to additional noise contributions due to correlations with the lensing potentials (22ll . 



4 



XX' 




(\uh) 






(ll,l2) 




TT 


(L 


■li) 




(L 


•I2) 




TE 


cos 2((^i| - 


Vi,)(L 


li) 


(L 


•I2) 




TB 


sin2((/i,, - 




1|) 







EE 


COs2((/3|| - 


Vi,)(L 


li) 


cos 2{ip\^ - 




I2) 


EB 


^m2{(p\, - 


'^l2)(L • 


li) 


sin2((^i, - 


</'b)(L 


I2) 


BB 


cos 2{^p^^ - 


W2)(L 


li) 


cos 2(^i| - 




I2) 



TABLE I: Window functions which appear in Eq. i ll lb and Eq. l ll2b ; here L = li + I2 



is the detector noise variance per steradian area for temperature (X - T) or polarization (X = E or B), and cr^ - OfKhm I V8 In 2 
is the effective beamwidth of the instrument calculated from its full-width half-maximum resolution 6f^yi,„. We will assume fully 
polarized detector, for which Iwj -We - Wb- 

Quadratic combinations of CMB fields can be used as estimators of the lensing potential field and hence the intervening 
projected mass between us and the last scattering surface. Furthermore, a CMB-field-squared map appropriately filtered in 
Fourier space can serve as an optimal estimator for the deflection field. Optimal filters for quadratic estimators have been 
designed lITsIl 

flfxx'(L) = ^^f^ r P^^X\h)X'\h)Fxx'i\uh) , (8) 
L J {2nY 

where X and X' can be T, E, and B. The normalization A^x' is chosen such that (c/xx'(L))cmb = d{V,) = L<p 



Axx'(L) = 



7T^fxx'(luh)Fxx'(h,h) 



(9) 



where 



rx'x'irxxtf a 1 \ rxx't^xx'if a 1 \ 

t-;, t-;, JXX'{h,h) - JXX'ihAl) 

FYY,(h h) = — flO) 

'AAV'l^'Z; (-^XXl(-^X'X't(--X'X't(-^XXl _ (C^XX'lc^XX'M ' ^ ' 

h h h h h ' 

where for XX' = TT, EE, BB, and TE, 

fxx'{\uh) = Cf 'Wxz'(li,l2) + Cf ^Wxx'{\uh\ (11) 

and for X = [T,E], X' = B, 

fxx'ihM) = Cf 'Wxx'ihM) + Cf^""' -Wxx'ihM). (12) 

The window functions Wxx' are given in Table U In Fig. [1] we show the input power spectrum of the lensed CMB fields, the 
reconstructed deflection field, and the corresponding Gaussian noise as obtained by using the estimators defined in Eq. ([8]). 
Note that the average ( )cmb denotes an ensemble average restricted only to different Gaussian realizations of the primordial 
CMB and instrument noise but assuming a fixed realization of the large-scale structure. The unmarked average, ( ), means 
the average over the primordial CMB field and the large-scale structure realizations as defined in Eq. For the purposes of 
estimating the large-scale structure in the real observable universe, it is essential to ensure that the estimators after appropriate 
averaging over realizations are truly unbiased for a typical realization of the primordial Gaussian CMB field. As we will see in 
the next section that in the presence of non-zero systematic contamination, the estimators are biased. 



III. INSTRUMENTAL SYSTEMATICS EFFECT ON THE LENSING POTENTIAL POWER SPECTRUM 

In this section, we use the quadratic estimators to derive the effects of instrumental systematics on the reconstructed lensing 
potential power spectrum. We show that instrumental systematics can introduce non-Gaussian correlations of CMB fields. The 
instrumental systematics-induced CMB trispectrum gives an extra contamination to the reconstructed deflection angle power 
spectrum. In subsection IIII A I we consider the simple case of the TT estimator to explain how systematics contaminate the 
lensing reconstruction process. In subsection llllBI we take the EB estimator as an example to show how instrumental systematics 
in CMB polarization measurements affect the lensing reconstruction. In order to compare the experimental requirements for 
primordial B-mode detection and lensing potential reconstruction, we also calculate the effects of instrumental systematics on 
B-mode detection for a given inflationary energy scale. 
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FIG. 1: Left panel: CMB power spectrum for the TT, TE, EE and lensed BB power spectrum. The lower (upper) dashed black line shows the 
temperature (polarization) noise for Expl. The lower (upper) dotted black line shows the temperature (polarization) noise for CMBPol-like 
reference experiment. Center panel: Gaussian noise for different quadratic estimators for the Expl. Right panel: Same as central panel but 
for CMBPol-like reference experiment. Note that the EB estimator has the lowest Gaussian noise, and may thus be considered as the best 
estimator among all the possible quadratic estimators. 



A. A simple example of temperature systematics 

We introduce the calibration parameter (gain fluctuation of receivers) fl(n) for temperature measurement, defined as: 

f"'"'{n) = [1 +a(n)]f' (ft). 



(13) 



If we assume that there is no correlation between the lensing potential 4> and the instrumental systematic fl(n), the power spec- 
trum of the lensed CMB temperature with systematics correction term can be found as (more details are given in AppendixlAll: 



J (2^0^ " ^^'-'M ' J (2^^l>~'.l^'. f^'-*'^-*'] (2^ 



(14) 



This result is given to linear order in the lensing-potential power spectrum Cy' and the gain fluctuation systematics power 
spectrum C"" . The last term represents the bias introduced by the calibration systematics. In the absence of systematic eff'ects, 
it is easy to prove that the deflection angle estimator is (iirr(L))cMB = L(p{L,) = dTiii^), as desired. But once the contribution 
from the calibration parameter fl(ft) is considered, one finds: 



<^f7T(L)>CMB = drril-) + 



AtAL) 



(15) 



i.e. in the presence of a(L) the estimator drAi^) is a biased estimator for the deflection field in Fourier space. Consequently, 
the deflection angle power spectrum (t/77-(L)t/7-7(L)) would be biased due to systematic contamination, and is given by (see 
Appendix lAl for details) 
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+i27Tfcl'''Cj'^'[6u{li' + IMh' + h) + 6D{h' + 1i)5d(1i' + h)] 
+K^,,-|/rr(li,li')/rr(l2,l2') + C^,%,|/rr(li, hO/rrCh, 1i')]5d(L + L') 

+ [CTfaa(lul2)faa(h', h') + C;';^^^,^U(h^l')faa(h, h') 
+qi:V,,'|/««(llJ2')/««(l2,ll')]5D(L + L')j 

(InfdDiL + L') Cf + Nfljj.{L) + N^^].^^(L) + N^^l^^(L) + . 



(16) 



where we define /fla(li,l2) = C^^ + CjJ. In the last line, the first term in the square bracket is the deflection angle power 

spectrum The second term is the so called Gaussian noise Nfjj-j{L) which gives the dominant noise contribution to the 

variance of the deflection power spectrum. The third term N^jjj(L) is the the leading order non-Gaussian noise which is 

first order in Cf^ and gives correction to the dominant Gaussian noise Nfjj.j{L). The forth term, Njjjj{L), is the leading 

order instrumental systematic contribution to the variance and is first order in C^'^. The Gaussian noise Njjjj.{L) and the first 

order non-Gaussian noise Njjjj{L) have been previously calculated in ifTsll and ll2lll . respectively. The systematic noise term 

N'"^jj.j{L) is a new contribution to lensing power spectrum. In principle, one should include noise terms which are higher order 

in C^'^ and C^'^ , however since both of them are small, we truncate at the first order and expect that higher order contributions 
are much smaller 

The quadratic estimator given in Eq. ^ is optimized in the presence of Njjjj(L), and assuming no contribution from the 
first order non-Gaussian noise N^jjj(L) and instrumental systematic noise Njjjj{L). Hence the estimator is optimal^ as long 
as A^^'^ jjiL) <K Nfj jj{L), and Njjjj{L) <si Nfj j-j(L). It has been shown that the non-Gaussian noise A^^'^ jjiL) is about one 



order of magnitude smaller than the Gaussian noise contribution lulh . In llllBI we calculate the systematic noise term (L) 
contribution to the lensing potential reconstruction for the EB estimator. 



B. General analysis on polarization systematics 

We parametrize the fields of instrumental systematics for CMB polarization measurements following ll32ll . The polarization 
contaminations fall into two categories, one associated with the detector system which distorts the polarization state of the 
incoming polarized signal (Type I hereafter), and another associated with distortion of the CMB signal due to the beam anisotropy 
(Type II hereafter). This parametrization can be generalized to different polarimeters. The instrumental response to incoming 
CMB radiation is usually described by the Jones transfer matrix. Bias induced in the matrix determination will mix the Stokes 
parameters determined from it. To first order, the effect of Type I systematics on the Stokes parameters can be written as ll32ll 



61Q ± iU]{h) = [fl ± i2co]mQ ± iU](h) + [/i ± //2](n)[e + iU]{h) + [y, ± iy.mmA). (17) 

a is a scalar field which describes the miscalibration of the polarization measurements (recall that in last subsection, we used 
a to denote the miscalibration of temperature measurements), oj is also a scalar field that describes the rotation angle of the 
instrument, (/i ± 1/2) are spin ±4 fields that describe the coupling between two spin states (spin-flip), and (yi + /ya) are spin ±2 
fields that describe monopole leakage from the temperature to polarization. 



^ If the systematic contributions are comparable to the Gaussian noise contribution then we need to design new optimal lensing reconstruction estimators to 
take into account the instrumental systematics effect. 
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Similar to the Type I systematics, the effect of Type II systematics on the Stokes parameters can be written as ||32|] 

6[Q ± iU]ih; cr) = crp(n) ■ V[2 + /t/](n; cr) + o-[di + id2]{n)[di + idoJTih; cr) + cr-^(n)[5i + id2fT{a; cr) (18) 

the systematic fields are smoothed over the average beam cr of the experiment. Therefore the type II systematic fields are 
sensitive to the imperfection of the beam on the scale cr. {pi + ip2) are spin +1 fields that describe pointing errors, (di + id2) are 
also spin ± 1 fields that describe dipole leakage from temperature to polarization, and qisa scalar field that describes quadrupole 
leakage iH. 

As a simple model, we will assume that the contamination fields, as defined in (fTTT i and ( fTSl l. are statistically isotropic and 
Gaussian (although some of the systematics fields need not be so), thus their statistical properties can be fully described by their 
power spectra, 

<5 (1)5(1')) = {Infdil + l')Cf^ , (19) 

where S stands for any of the 1 1 systematic fields. The systematic fields can be modeled with the power spectra of the form 

Cf^ =C„exp(-/(/+l)a2), (20) 

i.e. white noise above certain coherence scale as , which is a key quantity to affect the level of contamination of each systematics 
effects. The normalization factor Co can be determined by 



(21) 

where As characterizes the nns of the contamination field S . 

The instrumental systematics induce distortions on the CMB fileds. The contaminations to the BB and EE power spectra due 
to different measurement systematics take the form 

= Z / +g/ Uc^f;,C-(c^)[ly^(l,l')]^ (22) 

= Z / (02C^-^.'lC(-)[^l(l'l')]^ + g / (glC^-^.:|Cf[(cr)[lV|(U')]^ . (23) 

The explicit forms of W|(li,l2) and W^(li,l2) are given in Table HH which are the window functions of each systematic S 
for B-mode and E-mode harmonics, respectively. The summations of the first term on the RHS of Eq. (l22l i and Eq. ( |23l ) run 
over calibration a, rotation oj, spin flip /„ and fh, and pointing error ja and yi,. The summations of the second term run over the 
rest of the systematics parameters which describe the temperature leakage given in Table |ll] Cf^{cr) and Cj^{cr) are the beam 
smoothed temperature and E-mode polarization power spectra. 

Cf^(cr) = C,™ exp(-Z(Z + 1)0-), Cj^icr) = Cf^ exp(-Z(Z + l)cr). (24) 

We use Eq. ( l22l i and Eq. ( |23l ) to calculate the systematic requirements for B-mode detection in order to compare the require- 
ments for lensing reconstruction. We show the results in TablesHIJandllVI 

Now we move on to calculate the systematic contamination on the lensing power spectrum. The polarization fields are 
essentially uncorrected with the lensing potential, so if we do not consider secondary effects, the «-point functions with n odd 
are zero. The next non-zero order is the trispectrum. The calculation for the connected part of the trispectrum of polarization 
in the presence of instrumental systematics is similar to the connected temperature trispectrum presented in the last section. 
Here we give the results for the trispectrum related to lensing reconstruction using EB estimator, and refer the readers to the 
AppendixIClfor the explicit calculation for other quadratic estimators. At leading order we have 



d-l 

(2^ 



exp(-l(l + l)a'^) 
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Type of S 
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- sin[2((^;| - (pi) 


Dipole da 
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TABLE II: Window functions for all the 1 1 systematic parameters. First column indicates the type of systematic parameters in consideration. 
Second and third columns show window functions for systematics induced B-mode Wg(Ii,l2). and for E-mode W|(li,l2) respectively. These 
window functions are needed to calculate systematic contamination on primordial gravitational wave detection or the deflection angle power 
spectrum reconstruction. We note that I, = ZJi, I2 = L - li, and I2 = l^ii- 



(£(li)°''^B(l2)*^£(l'i)°''^B(l^)°''^)^ = (271)2^0(11 + I2 + I'l + I2) 



f^EE f^EE 



qf^i.iW^BCh, -11)1^^(12, -I'l) + C+i;|W^B(i2, -i'i)w^b(12, -ii) 

P-dislortion P— distortion 



s s 

T— leakage 

+ crcr[ Z c^f,,|W|(i2,-ii)w^(i;,-i;) + 



s 

T -leakage 



s 



(25) 



where we defined the lensing B-mode window function Wb(1, 1') = 1' ■ (1 - l')sm2(if\ - (p\>) and W|(li, -I'j) is the systematics 
window function for any of the 1 1 systematics parameters. The formula for each of the systematic window functions Wg (1, 1') 
can be found in Table HJ 

Different trispectrum can be constructed from combinations of the temperature and polarization fields. We discuss other cases 
in the AppendixICjin order to calculate EE, TE and TB estimators including contributions from systematics contamination. The 
formulas shown here are readily generalized to the full sky. For a discussion of the spherical generahzation of the polarization 
trispectra, see ll47ll . We have shown that E and B-modes are mixed not only by weak lensing, but also by instrumental sys- 
tematics. Even if there is no lensing induced correlation, certain kinds of systematics can give a non-zero contributions to the 
trispectrum. Now we move on to construct quadratic lensing estimators from E/B polarization modes, and we quantitatively 
show how instrumental systematics contamination affects the lensing reconstruction process. Again we take the EB estimator 
as an example, and leave the discussions of other estimators to the Appendix ICl The variance of the deflection angle power 
spectrum including systematic effects can be written as 



■ '^£b(L'))cmb), „„) 



Aeb(L)Aeb(L') 



~j (^2FEBih,h)FEBih,h) 



LSS/SYS 



L U 
x{E(li)°^'B(l2)°^'E(l[)°^'B(l'^r^') 



.(0) 



(1) 



j{S) 



(26) 



where L = li -H I2, and C'''^(L) is the deflection angle power spectrum. The terms Nfl ^g(L), N^^^ ^^(L), and N^^g ^^^(L) are the 
Gaussian noise, first order non-Gaussian noise, and the first order systematics noise. The Gaussian noise contribution comes 
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FIG. 2: Upper panels: Contamination from the calibration systematics a to the deflection angle power spectrum using the EB estimator. The 
rms fluctuation, Ag is assumed to be 10%. The left and right panel are for Expl and reference experiment respectively. In both the panels, 
the solid black, dashed blue, and dot red curves show deflection angle power spectrum C''''(L), Gaussian noise yV""(L), and the first order 
non-Gaussian noise ^"(L). The remaining curves show the absolute value of the systematic bias for various choices of coherence length a^, 
starting from as = 10' to a, = 120'. Lower panels: Same as the upper panels but for rotation systematics a>. 



from the disconnected part of the four-point function, while both the first-order non-Gaussian noise A^^'^ ^nd systematics 

we have used filters Feb which are given 



noise Np^„JL) contribution comes from the connected part. In deriving Eq. 



in Eq. (ITOl i. and the trispectrum which is given in Eq. (|25T l. The ellipses stands for terms beyond first order in the systematics 
power spectra or the deflection angle power spectrum. We note that the Gaussian noise term also includes systematic effects 
implicitly since instrumental systematics bias the measured power spectrum as we have shown. The Gaussian and non-Gaussian 
noise terms in the absence of systematic contribution have been previously reported in ifTsIl and ll2lll . The systematics noise 



(S ) 

term N^^^giL) is new, for which the explicit form is given as 



Aeb(l) Aeb{l') r d% r d% 



{L') r d% r d% , , 

■— ■ T -^FEBihMFEBih Ai ) 



P-distortion 



P-distorlion 



f-iEE^EE 



-I- cf^cF 

'1 I, 



2 C^f.,,|H^^(l2,-li)H^^(l^,-i;) + J] Cl^,v^iK(h,-l\)W',{l'„-h) 



s 

T -leakage 



T -leakage 



(27) 



Eq. dZTb is our main result for the systematics contamination on lensing power reconstruction using the EB estimator We will 
use this equation to numerically compute the systematic-induced bias for the 1 1 systematic parameters. The results are shown 
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FIG. 3: Upper panels: Contamination from the spin-flip systematics /„ to the deflection angle power spectrum using the EB estimator. The rms 
fluctuation, As is assumed to be 10%. The left and right panel are for Expl and reference experiment respectively. In both the panels, the solid 
black, dashed blue, and dot red curves show deflection angle power spectrum C''''(L), Gaussian noise ^""(L), and the first order non-Gaussian 
noise A'*"(L). The remaining curves show the absolute value of the systematic bias for various choices of coherence length o-,, starting from 
Q-j = 10' to ffj = 120'. Lower panels: Same as the upper panels but for spin-flip systematics ft- 



in Figures (O to (|7|i and are discussed in next section^. 



IV. RESULTS 



Figures|2]to|7]and TablesHIIl&IIVIsummarize our main findings. We have focused on the systematics-induced bias for the EB 
estimator (results for the other estimators are provided in the Appendix ICTi because it has the highest Gaussian signal-to-noise 
ratio for reconstructing the projected matter power spectrum (See Fig.[T]or |15il ). Figures|2]|2]show the contamination introduced 
by different systematic effects (the term N^^^^giL) in Eq. ( |27] i) in the deflection angle power spectrum reconstruction. We have 
assumed the rms fluctuation of the systematics fields to be 10%, and varied coherence length starting from minimum a, - 10' 
to the maximum as - 120'. For comparison, we also show the level of Gaussian noise and first order non-Gaussian noise of the 
EB estimator, the terms N^^^giL) and N^^^^giL) respectively in Eq. ( [27] ). Systematics-induced bias generally increases with 



^ As a note, for each instrumental systematics parameter S . for both the polarization distortion and temperature leakage, there are two terms which contribute 
to the final results. The terms proportional C^^^^| | is generally much smaller than the terms proportional to but for non-Gaussian noise M'' l2lll . the 

terms proportional to C?* . , are about an order of magnitude larger than other two terms proportional to C^^ , and C^* , 
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FIG. 4: Upper panels: Contamination from the pointing systematics p„ to the deflection angle power spectrum using the EB estimator. The 
rms fluctuation, As is assumed to be 10%. The left and right panel are for Expl and reference experiment respectively. In both the panels, 
the solid black, dashed blue, and dot red curves show deflection angle power spectrum C''''(L), Gaussian noise N^^^L), and the first order 
non-Gaussian noise N^^^(L). The remaining curves show the absolute value of the systematic bias for various choices of coherence length a„ 
starting from o-j = 10' to o-, = 120'. Lower panels: Same as upper panels but for the pointing systematics pi,. 



the increase in coherence length. For all the coherence lengths, the bias N^g^g(L) is fairly constant on large scales, similar to 
the feature of the Gaussian noise and the first order non-Gaussian noise of the estimator, however on small scales [i > 200), the 
systematic contamination are not constant, and for some systematics, they oscillate between positive and negative values. The 
fluctuating features are due to the oscillation of Gaussian noise at high i and the CMB power spectra included in the systematic 
calculations, smoothed and filtered by certain window functions. The negative systematic contributions at certain multiple range 
for some parameters are caused by combining the eff'ects from window functions, which sometimes give negative, and TE cross 
correlation which enters the calculation of monopole leakage, dipole leakage, and quadrupole leakage from temperature field. 

In Table Hn] we calculate the maximum rms amplitude As required for each systematic parameters to keep its own contami- 
nation lower than the deflection angle power spectrum. Lensing extraction requires the control of systematic rms fluctuations 
at levels depending on the type of systematics. Among all the 1 1 systematic parameters, rotation w, and monopole leakage 
ya&yb place the most stringent requirements, while pointing error pa&Pb, quadrupole leakage q, and calibration a are among 
the least demanding. Also with respect to sensitivity to the coherence length, some systematic parameters such as rotation w and 
monopole leakage ja&7h are very sensitive (as shown in Figure ^ and (|5j), while some are not very sensitive. 

For comparison we also show the systematics requirements for primordial B-modes detection. For B-modes we consider two 
inflationary energy scales £ = 1.0 x 10'^ GeV and E = 3.0 x 10'^ GeV. In addition we consider two choices of coherence length 
Us - 10' and o-j = 120', and two experimental setups Expl and reference experiment. Since the systematic requirements from 
B-modes detection is sensitive to the inflationary scale in question, it is useful to define an energy scale Ea-ii below which lensing 
reconstruction is safe and above which lensing sets the systematic requirements. We define the critical energy of a given 
systematic parameters as the inflationary energy for which the systematic requirements to extract the B-modes signal is equal to 
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FIG. 5: Upper panels: Contamination from the dipole leakage da to the deflection angle power spectrum using the EB estimator. The rms 
fluctuation, As is assumed to be 10%. The left and right panel are for Expl and reference experiment respectively. In both the panels, the solid 
black, dashed blue, and dot red curves show deflection angle power spectrum C^'CL), Gaussian noise A'""(L), and the first order non-Gaussian 
noise A'*"(L). The remaining curves show the absolute value of the systematic bias for various choices of coherence length o-,, starting from 
Os = 10' to Q-j = 120'. Lower panels: Same as upper panels but for the dipole leakage di,. 



the requirement to extract the lensing signal. Hence for a given systematic parameter, if the requirement for detecting B-modes 
with energy scale E < Ecrn. have been met, then lensing extraction requirements are already met. Similarly if a given systematic 
parameter sets the detectable B-modes with energy scale E > Ecru., then lensing extraction requirements provide more stringent 
constraints on that systematic parameter. We show the critical energy Ecnt- in Table IIVI for all the 1 1 systematic parameters for 
two choices of coherence length (or, = 10' and Os - 120') and for both the experimental setups (Expl and reference). 

For both the Expl and reference experiment, lensing reconstruction is safe once the experimental requirements to detect B- 
modes with inflationary scale £, - 1.0 x lO'*" GeV are met. The requirements from lensing extraction are about 1-2 orders of 
magnitude less stringent than the requirements to measure the primordial B-modes with inflationary energy scale of 1.0 x 10'^ 
GeV. This means that once the experiment satisfies the requirements of systematic control to detect primordial B-modes, it is safe 
to use such an experiment to reconstruct the lensing potential power spectrum without extra effort on improving the instrumental 
systematic control. On the other hand the requirements for lensing reconstruction are comparable or even more stringent for 
some systematic parameters than the requirements to detect primordial B-modes with inflationary scale £, = 3.0 x 10'^ GeV. 



V. CONCLUSION 



We illustrate the eff'ects of instrumental systematics on the reconstruction of the projected matter power spectrum from CMB 
gravitational lensing. We consider seven types of effects which are related to known instrumental systematics: calibration, 
rotation, pointing, spin-flip, monopole leakage, dipole leakage and quadrupole leakage. These effects can be parametrized by 1 1 
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FIG. 6: Upper panels: Contamination from the monopole leakage ja to the deflection angle power spectrum using the EB estimator. The rms 
fluctuation, As is assumed to be 10%. The left and right panel are for Expl and reference experiment respectively. In both the panels, the solid 
black, dashed blue, and dot red curves show deflection angle power spectrum C^'CL), Gaussian noise A'""(L), and the first order non-Gaussian 
noise A'*"(L). The remaining curves show the absolute value of the systematic bias for various choices of coherence length or,, starting from 
Q-j = 10' to Q-j = 120'. Lower panels: Same as the upper panels but for the monopole leakage yi,. 



distortion fields. Each of these systematic effects can mimic the effective projected matter power spectrum and hence contaminate 
the lensing reconstruction. We assume a Gaussian distribution for each parameter. We have modeled the fluctuations in the 
instrumental contamination fields with a coherence length (o-j) and an rms amplitude (A,), as defined in Eq. ( |20] i and ^T\\ . 
respectively. Rotation systematics o) and monopole leakage ya&y/, are among the most dangerous ones. 

It is important to know how the systematic effects propagate to the lensing potential reconstruction with good precision, 
in order to reliably reconstruct the lensing potential. Without well understood and caUbrated instrumental systematics, the 
extraction of cosmological parameters from the statistical properties of lensing potential would also be biased. Systematics 
in the lensing reconstruction are also one of the key concerns for detecting primordial gravitational waves by removing the 
contribution properly from lensing-reduced B-modes. A faithful detection of the primordial gravitational wave signal largely 
depends on how clean the delensing process has been. 
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FIG. 7: Contamination from the quadrupole leakage q to the deflection angle power spectrum using the EB estimator. The rms fluctuation. As 
is assumed to be 10%. The left and right panel are for Expl and reference experiment respectively. In both the panels, the solid black, dashed 
blue, and dot red curves show deflection angle power spectrum C''''(L), Gaussian noise N^'^\L), and the first order non-Gaussian noise A'''*(L). 
The remaining curves show the absolute value of the systematic bias for various choices of coherence length a^, starting from = 10' to 
a, = 120'. 
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a, = 10' 
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a, = 10' a, = 120' 


a, = 10' 


a, = 120' 


a, = 10' 


a, = 120' 


a, = 10' 


a, = 120' 


a, = 10' 


a, = 120' 


Calibration a 


1.04 


0.55 


0.549 


0.468 


0.061 


0.052 


0.86 


1.56 


0.486 


0.468 


0.054 


0.052 


Rotation w 


0.30 


0.061 


0.27 


0.207 


0.030 


0.023 


0.24 


O.Il 


0.243 


0.198 


0.027 


0.022 


Pointing 


2.10 


1.57 


8.55 


6.12 


0.95 


0.68 


1.93 


5.55 


15.3 


12.24 


1.70 


1.36 


Pointing 


1.39 


2.76 


1.08 


6.57 


0.12 


0.73 


2.31 


9.74 


1.71 


13.14 


0.19 


1.46 


Flip/,, 


1.08 


0.17 


0.549 


0.441 


0.061 


0.049 


0.62 


0.31 


0.486 


0.432 


0.054 


0.048 


Flip .A 


0.61 


0.17 


0.531 


0.360 


0.059 


0.040 


0.58 


0.31 


0.477 


0.36 


0.053 


0.040 


Monopole y„ 


0.114 


0.024 


0.021 


0.0058 


0.0023 


0.00064 


0.13 


0.013 


0.0207 


0.0058 


0.0023 


0.00064 


Monopole yb 


0.114 


0.036 


0.014 


0.0034 


0.0016 


0.00038 


0.64 


0.12 


0.0144 


0.0034 


0.0016 


0.00038 


Dipole du 


0.82 


0.11 


0.085 


0.060 


0.0094 


0.0067 


0.97 


0.39 


0.153 


0.117 


0.017 


0.013 


Dipole db 


0.55 


0.092 


0.85 


0.063 


0.0094 


0.0070 


1.02 


0.33 


0.153 


0.126 


0.017 


0.014 


Quadrupole q 


1.58 


0.78 


0.162 


0.558 


0.018 


0.062 


2.38 


5.47 


0.495 


2.25 


0.055 


0.25 



TABLE III: Systematic contamination for lensing and B-mode detection for the two experimental setups (Expl and reference), for two choices 
of coherence length a, = 10' and a, = 120'. We scale the nns amplitude of the systematics field A, to the same level as the maximum of the 
signal (deflection angle power spectrum at f = 40 and B-mode signal for a fiducial energy scale at £ = 90). First column mentions the type 
of systematic parameter in consideration. Within the lensing columns we show the maximum A, (for as = 10' and a, = 120') for bias to not 
exceed the lensing power spectrum. For B-modes columns we show the maximum A, (for a, = 10' and as = 120') for bias to not exceed the 
B-mode signal for two choices of inflationary energy scale £ = 3.0 X 10'* GeV and £ = 1.0 x 10'* GeV respectively. 
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Type 
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Reference 


£,„■,. (GeV) 


£„-,Y.(GeV) 


a, = 10' 


a J = 120' 


= 10' 


ff, = 120' 


Calibration a 


4.13 


3.25 


3.99 


8.83 


Rotation vi' 


3.16 


1.63 


2.98 


2.24 


Pointing p„ 


1.49 


1.52 


1.07 


2.02 


Pointing pi, 


3.40 


1.94 


3.49 


2.58 


Flip,/;, 


4.21 


1.86 


3.34 


2.54 


Flip/z, 


3.22 


2.06 


3.31 


2.78 


Monopole y„ 


7.80 


5.00 


7.80 


4.33 


Monopole 74 


8.29 


9.73 


20.00 


17.77 


Dipole d„ 


9.34 


4.05 


7.55 


5.48 


Dipole dt, 


7.51 


3.57 


7.78 


5.00 


Quadmpole q 


9.40 


3.48 


6.02 


6.42 



TABLE IV: We show the critical inflation energy scale E„ii. for each systematic parameter, for two choices for coherence length (o-, = 10' and 
Q-j = 120') and for two experimental setups (Expl and reference). First column mentions the type of systematic parameter in consideration. 
Ecrii of a given systematic parameter is the inflationary energy for which the systematic contamination (from that parameter) is equal to both the 
maximum B-modes signal (corresponding to ^ = 90) and the lensing signal (corresponds to £ = 40). Hence for a given systematic parameter 
if the requirement for detecting B-modes with energy scale E < E„ii. have been met, then lensing extraction requirement are already met. 
Similarly if a given systematic parameter sets the detectable B-modes with energy scale E > , then lensing extraction requirement provide 
more stringent constraint on that systematic parameter. 
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APPENDIX A: INSTRUMENTAL SYSTEMATICS FOR THE TT ESTIMATOR 



In this appendix, we show detailed calculation of instrumental systematic contamination on the TT estimator. The calcu- 
lation is similar to the systematic contamination analysis for other estimators involving the polarization fields (discussed in 
Appendix IbTi. 

We first expand the observed temperature T"'^%n) field, given by in Eq. ( fT3T l as 

= [1 +a(n)]r[n + V0(n)] 

^ Tin) + V„cf,{nW"T{n) + iv„,^(n)Vi0(n)V°V''r(n) + . . . 

+a(n) X [Tin) + V„^(n)V°r(n) + ^V„cf>{hWb4>{nW"^'T{n) + ...], (Al) 

where fl(n) is the calibration parameter. Taking the Fourier transform of the lensed map with systematics 7'"'"(n) under the 
flat-sky approximation, 

f°^'(l) = J £/fi[l -Hfl(n)]f(n)e"'' * 

= m- f 4^r(l')L(l,l'), (A2) 

where 

L(l,l') = 0(1-1')[(1-1')-1'] + ^ r ^0(1") X 0(1-1' -l")(r-l') [(!" + !' + ■■■ 

-fl(i-i')+ r -^fl(r)x[i'-(i-i'-r)]0(i-i'-i") 

A [ Fr^«('"') ^ "^(1") ^ "^(l - 1' - 1" - • >') + 1" + 1' - 1) • 1'] + . . . . (A3) 

2 J (2n)- J (IttY 

Using the above equations one can compute the CMB temperature power spectrum Cj^ in presence of systematic contamination, 
the result is given in Eq. (fT4l i. 
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Analysis of the deflection angle power spectrum Cf'^ involves trispectrum calculation. Below we show the details of the 
temperature trispectrum including the systematic contribution. 



\{C]J + C,^^") (27r)25D(L) + fl(L) [C]J + C,7) + 0(L)/rr(li,l2) - J ^fl(li')(L - li')0(L - liOhCf 

- J |^fl(l2')(L - l2')0(L - l2')liCf + i J |^0(li')0(L - li'){ci;^(li ■ li')[li ■ (L - li')] + Cl\h ■ h')[h ■ (L - li')]} 
+ r -li>(l2 + li')Cr- r ;^C,r0(li-li')0(l2 + li')[li'-(li-li')]x[li' -(12+11 ')] 

- J ^C,r,^(l2-l2')fl(ll+l2')[l2'-(l2-l2')]- J ^C,T^0(ll-llXl2 + ll')[ll'-(ll-ll')]| 

(C,r + C,T^") {Infs^d.) + fl(L) (C,T^ + CjJ) + 0(L)/rr(li', I2') - J (^«(1)(L " D-^d^ - Dh'Qf 

- r -^fl(l)(L-l)0(L-l)li'C,T^ + i r ^0(l)<^(L-l){c,r(li'-l)[li'-(L-l)]+C,r(l2'-l)[l2'-(L-l)]} 

J (/jr)" 1 2 J (2;r)^ ^1 2 ) 

+ r 7^«(li'-lMl2'+l)Cr- r -^Cf^0(li'-l)0(l2'+l)[l-(li'-l)]x[l-(l2'+l)] 

J (27r)^ J (2;r)- 

^Cj^4>(h' - + 1)[1 ■ (I2' - 1)] - J (^Cr<^(li' - l)«(l2' + • di' - m+Penn. (A4) 

In practice, we do not know the large-scale structure between us and the last-scattering surface, we must average over diff'erent 
realizations of the large-scale structure (denoted by ( )lss) and also over different realizations of the systematic fields (denoted 
by ( )sYs) to obtain the observed deflection angle power spectrum. The deflection field power spectrum can be estimated by 
the TT estimator as 

((<Jrr(L) ■ flfrr(L'))cMB)^^j, J^YS = ^"^^^^^^ ' '^rr(L')) 

a J ,^ Att{L) Att{L') 

x/l^/ ||((<7'°Mi)r°M2)rM:')rM/)>cMB),,,),,, 

xFTT(h,h)FTTih',h')- (A5) 

Only the connected part of the trispectrum comes into the deflection angle power spectrum C^''. Keeping the first order of Cj*"* 
and C°", the results can be simplified as 

{r'\h)T'''\l2)T'''\h')T"'\l2')}cMB,conn. = ( j ^T{h')4>(h - h') [{h - h') ■ h'] 

X J ^r(i2')<^(i2-i2')[(i2-i2')-i2']r(ii')r(i2')) 

+( J ^r(ii')fl(ii-ii') J ^r(i2')fl(i2-i2')r(ii')r(i2')) 
= cj,cl{<p(h + h'mi + 12')) [di + ii') ■ ii'] [(I2 + 12') ■ I2'] 

+Cj:cl{a{h +li')a{l2 + I2')) 

+cici{<p(h + h'mi + ii')) [di + 12') ■ 12'] [(I2 + ii') ■ ii'] 

+Cj:cJ^{a(li + hXh + li ')) + Perm. , (A6) 

The terms that are linear in the lensing potential 0(1) and the temperature systematics fl(l) vanish. Further averaging over 
the lensing field and the systematics field, one can get Eq. ( fTSI ) which gives the measured lensing induced CMB temperature 
trispectrum in the presence of instrumental systematic contaminations. 
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APPENDIX B: POLARIZATION SYSTEM ATICS CONTAMINATION 

In this Appendix we calculate the window functions (given in Table ^ for each instramental systematic. We show how 
the systematic contamination comes in the analysis inevitably. Again the calculation is done for small deflection angles (weak 
lensing). 

We will first do the calculation for the polarization transfer systematics. The Taylor expansion of two Stokes polarization 
parameters including the systematic parameters can be written as: 

[Q ± iUr'-'in) = [l+a± ilumm ± iU][n + V,^(n)] + [/i ± ifiKmQ + iU][h + V<p(n)] + [yi ± in](n)m + V0(n)] 
= [g ± iU]{n) + V„0(n)V"[e ± iU](n) + iv„0(n)Vi0(n)V«V''[e ± iU]{n)] + ... 

+ [a ± i2co]m [Q ± iUm + V„<^(n)V"[e ± iU]{n) + iv„0(n)Vi0(n)V"V''[e ± /f/](n) + . . .] 

± i72](n)[ [Q + iU](n) + V,,^(n)Vne + iUm + ^V„0(n)V,<^(n)V"V''[e + iUm + ...]] 

+[ri ± iyojmTih) + V„0(n)V"r(n) + iv,0(n)Vfc0(n)V"V*r(n) + ...]. (Bl) 

Performing the harmonic transformation we get 

[E(l) + iBil)]"''' = J t/n [g + /f/]°*'(n)e'-''^'e"'' * 

= E(l)±iB(l)- f ^m')±iB(V)]Lp(l,V)- f -^r(l')Lr(l,l'), (B2) 

where 

Lp{\,\') = _ 1') [(1 _ r ^e±2,(^,,-<„)^(,") X 0(1 _ 1' _ 1") (1" . 1') [(1" + 1' _ 1) . 1'] + . . . 

2 J {2n)- 

e±^i(fv-f^)^a(\") + ;2w(l")] X [1' ■ (1 - 1' - 1")] 0(1 - 1' - 1") . • . 

_g±2,(2^,_,,-.„,-.„)[yj(i _ 1') + ,^2(1 _ 1')] + r lIle^2/(2^,„-.„,-.„)[yj(i") + ^^^(l")] X [1' ■ (1 - 1' - 1")] 0(1 - 1' - 1") . . . 

J {27lY 

(B3) 

and 



2yi 

±2'Xw-i'-w)r^,, n _ V\ + ,^/„n _ I'm j. I Jf ^±2/(^i-i'-w)r^,. n"^ 

(B4) 



Lr(l, 1') ^ -e^2,(^.-.'-^')[yi(l - 1') ± i7i{\ - 1')] + r ;^e^"(^'-"-^''[ri(l") ± /yzd")] x [1' ■ (1 - 1' - 1")] 0(1 - 1' - 1") 



Similarly for the local coupling systematics, the Taylor expansion of two Stokes polarization parameters including the sys- 
tematic parameters can be written as 

[Q ± iU]''''\n) = (1 + 6)[Q ± iUW + V0(n); cr) = [Q + /t/](n + V0(n); cr) + crp(n) ■ V[e + /;7](n + V0(n); cr) 
+(r[di ± id2\{rL)[di ± idiVin + V0(n); cr) + crq{rL)[dx ± idifTin + V0(n); cr) 

^ [Q± /t/](n; a-) + V„0(n)Vne ± /f/](fl; cr) + iv„0(n)Vi0(n)V«V*[e ± jf/](n; cr)] + . . . 

+ap(n) ■ V[ [2 + /f/](n; cr) + V„0(n; (r)V''[e + /f/](n; cr) 

+ iv,0(n; cr)Vi0(n; cr)V^V*[e ± jf/](n; cr) + . . .] 

+o-[£/i + /£/2](n)[5i + !52][r(n; cr) + V«0(n)V"r(n; cr) 

+ iv„0(n)V,,0(n)V"V''r(n;(r) + ...]. 

+cr2^(n)[5i + id2\\T{n- a) + V«0(n)V''r(n; cr) 

+ ^V,0(n)Vfo0(n)V'^V''r(fi; cr) + . . .] . (B5) 
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Again performing the harmonic transformation gives 

[E(l) + iB(l)]"'" = J dn[Q±iU]''''\n)e^^'^'e-'^" 

= £(1)±/B(1)- r -^[£(l')±/B(l')]Mp(l,l')- r 4^r(l')Mr(l,l'), (B6) 

where 

Mp(l, 1') = _ 1') r(i _ 1') . I'l + i r .^e±2,(w-w)^(i")0(i _ i' _ i") (j" . j') r(i" + i' _ i) . i'] + _ , _ 

+/crl' ■ [/?„(! - 1') + ipi,{l - l')]e±2,(;.,,-;.,)g±'V,-,' 

^''^ r ^ ± • (1 - 1") 1")] -^(1 - 1' - . . . , (B7) 

and 

J'-)]g±'(</'i'+i^i-i'-2¥'i) 

" r ^-T? ± '^^^^(1")] ■ [e^''^'-''-'" (1 - r - 1") + e^'-f' 1'] [1' • (1 - r - 1")] </i(l - 1' - e^^'"^' . . . 

+0-2^(1 -l')l''e='-'*^''"'^'' 



Mr(l,l') = !o-l'[c/«(l-l')±iWi(l-l')]e^ 
-icr 



r2 f — 

' j (2;r)- 



2?(l")[e^'''*'^'-''-'""'^''(l - 1' - 1")^ + e±2;(w-^,)i/2-| .(!_!'_ j")] _ j' _ j") _ _ (B8) 



Using equations (IB21 l- (IB8b . the window functions for each systematic fields (defined in Table can be obtained without any 
difficulty. 

APPENDIX C: SYSTEMATIC CONTAMINATION IN OTHER ESTIMATORS 

In order to keep the paper short and clean, we focused on the EB estimator, one among six possible quadratic estimators of 
CMB fields. In this appendix, we give the results for EE, TE, and TB estimators. 

1. EE estimator 

For CMBPol like experiment, for multiple / < 100, EE estimator has comparable noise for the C'j'' reconstruction (especially 
in case of noise dominated B-mode detection). If one has good signal to noise ratio on E-mode detection, EE estimator can be 
a better choice compared to EB estimator Here we calculate the systematic contamination for EE estimator. Again, we first 
calculate the connected trispectrum, which can be written in a compact form using the window functions defined in Table HIl 

{E{hf'°'m2f'''E(\'if'°'E{V^r'°')^ = {2nf6u{h +h+ \[ + Ij) x 

[Cf C,f (^dl + + l^))W£(l2, -h)WEil'2, -I'l) 

P-distorlion 

+ Yj Cl''Cf{s (li + h)S 'd'l + l^))H'|(l2, -li)W^|'(l2, -I'l) 

SS' 
T-!eakage 

+ Yj Cl^Cl^{Sih + l2)5'(l'i + 12)>W|(12, -li)W|'(l2, -I'l) 

SS' 

+Cf C,f (,^(li + l'2)<p(h + I'l)) WeHi, -1[)We(1'2, -li) 

dislorlion 

^ Cf (5(li + 1^)5 '(l2 + I'l)) W'^ih, -l[)Wl'(l'2, -li) 

SS' 
leakage 

+ 2 Cl^Cj:f{S(h + l2)5'(l2 + I'l)) ^^£(12, -l'i)W^f (l2, -li) + Perm. I (CI) 



P-distoriion 

+ 

SS' 
T -leakage 
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Further averaging over the lensing and the systematic field gives 
(£(li)*^£(l2)°''^£(i;)°''^£(i;)*^)^ = (27Tf6u{h +h+ I'l + 12) X 



P-distortion 



'^^121^^(12, -\i)We{\2, -I'l) + q%|W£(l2, -ll) 

n P-distorlion 



+ '-1: 



T -leakage 



T-leakage 

s 



+ Perm , 



(C2) 



where we define the lensing E-mode window function as We(\, 1') = -1' ■ (1 - 1') cos 2{ip\i - The systematic window functions 
W|(li, -I'j) for any of those 11 systematics parameters are defined in Table HIl Hence we obtain the CMB E-mode polarization 
trispectrum due to gravitational lensing in the presence of systematic contamination. We define 



fp (ll , I2) = Cf , I2) + Cf; Wl (I2, ll), 



EE tjltS I 



and 



/|(ll,l2) = 1^1(11,12) + CfW|(l2,ll) . 



(C3) 
(C4) 



Thus Eq. ( |C2| i could be written as 
(£(li)°'''^£(l2)°'"^£(l'i)*-'£(i;)°'''^)^ = 



Cl^,^xjEE{\uh)fEE(h' ,h') + C*f^,,/££(ll,ll')/££(l2,l2') + Cj,%,|/££(ll , hO/ffiCh, h ') 



P-distortion 



Z 

55' 

T-leakage 

z 

SS' 



q;,ti/i£(u')/i£(u') + q;,,'|/i£(ii,ii')/i£(i2,i2') + ci^,,,fEEih,h')fEEih,h') 



q^,|/i£ai')/i£(u') + q,;vi/i£(iiJi')/i£(i2,i2') + ci'^,^jk(\uh')fEEihM) 



(C5) 



where the permutations now contain 5 additional terms with the replacement of (13,14) pair by other combination of pairs (li , I2), 
(li,l3), (li,l4), (12,13), (12,14)- One can simply Plug in the trispectrum given above to the deflection angle power spectrum 
estimation (the final result can be written as a sum of deflection angle power spectrum and noise terms) 



x(£(li)°''^£(l2)°'''£(l'i)°''^£(l^)°''^) 
= (27r)25D(L + L') 



(C6) 



where L = li + I2, and C'^'^{L) is the deflection angle power spectrum. The terms A^^^g^(L), A^^'^g^(L) and N^^^ ^^(L) are 
the Gaussian noise, first order non-Gaussian noise, and the first order systematics noise for the EE estimator respectively. The 
systematics noise term n''^^^^(L) can be written as 

l2 n\ n ^2i. r dH', 



N' 



L2 



r d%_ r £1 

J '(2^ J (2^ 



FEEih,h)FEEih' , h') 



P-distortion 



P-distortion 



r^EE r-^EE 



. f~^TE f~^TE 



Z C^f.,,|W^|(l2,-ll)W'£(l2>-l'l) + Y. C^'lil^£(l2'-l'l)^£(l2>-ll) 

5 5 
T-leakage T-leakage 

' Yj Cl^,i.jW',(h,'h)Wl(V„-l[) + Y C^f,,-|H'f(l2,-l'i)H'|(l2,-li) 



(C7) 



where the systematics window functions W| (li, I2) are given in Table HH 
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2. TB and TE estimator 

The instrumental systematic contamination calculation for the TB and TE estimator is similar to the calculation of EB estima- 
tor At leading order of Cj*^ and Cf^ ,we have 



(f(ll)°''^l(l2)°''^f(l'l)*^l(l2)°''^)^ = (iTtfSuih +h+l[ +I2) X 

jc,^^c,T^(0(ii +i2)0(i'i +r2))wx{h,-h)Wx{r2,-i\) 

P-distortion 

55' 
T-leakage 

+ J] Cl'^CmSdi + h)S'{l\ + 12))W^(12, -I'l) 

SS' 

+Cj^C]7{(l>{h + 12)<^(12 + I'l)) Wxih, -l'i)Wx(l^, -li) 

P-disiortion 

+ J] Cl^Cj/{S{h + 12)5 '(I2 + I'l)) w'^ih, -i\)w'^'(i'2, -ii) 



55' 

T-leakage 



+ ^ CfC,T^(5 (li + l^)5'(l2 + I'l)) W'^ih, -l'i)H^|'(l2, -li) 



= (2nf6u(h +h+ I'l + I2) X 



^TE^TE 



cll^^_^Wx{h, -h)Wx{\'2, -I'l) + qf^,,|Wx(i2, -i'i)Wx(i^, -ii) 

P-distortion 



P-distortion 

+ Z C"+l.l^x(l2,-ll)Wl(l^,-l'i) + C^f.lil^l(>2'->'l)^l(l2.-Il) 

s s 

T— leakage 

Z C^f.I.|W'l(l2,-ll)W'l(l2--l'l) + 



, y^TTy^TT 



T— leakage 

Z c^'m^l(»2,-i'i)H^^(i2,-ii) 

5 



(C8) 



where X - E for the estimator, and X - B for the estimator. We define the lensing window function ^^(1,1') = 
V ■ (1-1') sin2(i^i - if\>) and We{1, 1') = -1' ■ (1-1') cos 2{(p\ - ipy). Window W|(li, is the systematics window function for any 
of those 1 1 systematics parameters we defined in Table HIl One can simply Plug in the trispectrum given above to the deflection 
angle power spectrum estimation (the final result can be written as a sum of deflection angle power spectrum and noise terms). 



((<Jrx(L) ■ dTx(.V))cMB)^J 



Atx(L) AjxiL') 



LSS/SYS 



(L') r d\ r d% p 



L2 L 

x(f(ii)*^i(i2)°'''f(i'i)°''''i(i;)*') 

(2nf6D(L + L') 



(C9) 



where L = li + I2, and C'''^{L) is the deflection angle power spectrum. The terms Nf^j-^iL), Nj^^j^{L) and Nj^j-^iL) are the 
Gaussian noise, first order non-Gaussian noise, and the first order systematics noise for the TB/TE estimator respectively. The 
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systematics noise term A^' ' jx^^^ written as 



r d% r d% 

j (2n)2j (2n)2 



^Frxih , WFrxih ',h') 



P-dislortion 



P— distortion 



Z clU_\'^'xih,-h)wl{\',,-\\) + C^f.i;i^x(i2,-i'i)wl(i2,-ii) 

5 S 
—leakage T -leakage 

Z C^f...|W^l(l2,-ll)W^l(l2--l'l) + C^f„,|W^(l2,-i;)lV|(l^,-ll) 

S S 

where the systematics window functions W^(li,l2) are given in Table Ull 



(CIO) 



